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ABSTRACT 

We constrain blastwave parameters and the circumburst media of a subsample 
of BeppoSAX Gamma- Ray Bursts. For this sample we derive the values of the 
injected electron energy distribution index, p , and the density structure index of 
the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical 
and nIR afterglow data. The spectral fits have been done in count space and 
include the effects of metallicity, and are compared with the previously reported 
optical and X-ray temporal behaviour. Assuming the fireball model, we can find 
a mean value of p for the sample as a whole of 2.035. A statistical analysis of the 
distribution demonstrates that the p values in this sample are inconsistent with 
a single universal value for p at the 3cr level or greater. This approach provides 
us with a measured distribution of circumburst density structures rather than 
considering only the cases of k = 0 (homogeneous) and k — 2 (wind-like). We 
find five GRBs for .which k can be well constrained, and in four of these cases 
the circumburst medium is clearly wind-like. The fifth source has a value of 
0 < k < 1, consistent with a homogeneous circumburst medium. 

Subject headings: gamma rays: bursts 


1. Introduction 

Since the discovery of GRB afterglows in 1997 (Costa et al. 1997; Van Paradijs et al. 
1997; Frail et al. 1997), their relatively longer duration and broader wavelength coverage com- 
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pared with the prompt emission has made them the most accessible and arguably the most 
profitable area for observational GRB studies. Their non-thermal spectra and lightcurves 
can generally be well described by the fireball model (Cavallo & Rees 1978; Rees & Meszaros 
1992; Sari et al. 1998) which describes a decelerating, highly relativistic outflow, the so-called 
blastwave, interacting with a surrounding external medium. Adaptations have been made 
to accommodate a non-uniform external medium (Chevalier & Li 1999) and the fact that 
the outflow is collimated (Rhoads 1997, 1999; Sari et al. 1999). The evidence for collimation 
comes from energy considerations, and directly from observations of an achromatic break in 
afterglow lightcurves, the so-called jet-break. The jet-break is an indication of the time at 
which sideways spreading of the jet begins to become important, combined with the fact that 
the effect of relativistic beaming starts to be noticeable within the jet. For the long-duration 
GRBs, which most probably originate in the core-collapse of a massive star (Woosley 1993; 
MacFadyen et al. 2001; Woosley et al. 2002; Hjorth et al. 2003), a n(r) oc r” 2 density profile 
from a stellar wind would not be unexpected for the external medium. Its density can be 
probed through the temporal decays of the afterglow and, contrary to expectations given 
their stellar progenitors, long GRBs studied thus far do not all require a wind-like n(r ) oc r -2 
density profile in the local external medium (e.g. Panaitescu & Kumar 2002). 

The synchrotron nature of the blastwave produces spectra and lightcurves comprising 
of a set of power laws with characteristic slopes and frequencies. Accurate measurements of 
these observed quantities allow the underlying parameters of the blastwave to be determined 
(e.g. Wijers & Galama 1999; Panaitescu & Kumar 2001). Optical through X-ray spectra 
provide the opportunity to measure the index of the input power law energy distribution of 
electrons, p, potentially the index of the density profile of the circumburst medium, n(r) oc 
r _fc , and in some cases also the cooling break frequency v c — the frequency of electrons 
whose radiative cooling time is equal to the dynamical timescale of the blastwave. 

If the microphysics of all GRBs is assumed to be the same, the range of values of the 
input electron energy distribution indices should be narrow. However, the directly measur- 
able data that lead to the parameter p, such as the break frequencies and power-law slopes 
of the spectra and lightcurves, are dependent on various other factors, like the circumburst 
density, the fraction of energy in electrons (e e ) and magnetic fields (e s ), and simply the 
total amount of energy, making it more difficult to obtain a consistent value for p. Further, 
determination of these parameters so far works most successfully for the later-time afterglow 
considered here: before approximately 0.1 days the observed lightcurves and spectra are 
likely some mixture of the prompt emission (attributed to internally colliding shocks) and 
afterglow (the external shock), see e.g. Nousek et al. (2006) and O’Brien et al. (2006). These 
authors discuss the ‘canonical’ X-ray lightcurve for Swift bursts, which has three phases, the 
final phase beginning between 0.01 and 0.12 days after trigger (Nousek et al. 2006) and 
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showing the type of steady decay seen in the pr e-Swift era. The advantage of this is, though, 
that we can reliably restrict ourselves to the slow cooling regime in modeling the blastwave, 
where the injection frequency of the electrons is well below their cooling frequency. The 
measurements we perform in this study are made at times of approximately 1 day (see Table 
1, Starling et al. 2006, hereafter Paper I), so should not be affected by any prompt emission 
components. 

A potential problem is that occasionally values of p below 2 have been found (e.g. for 
GRB 010222, Masetti et al. 2001; Stanek et al. 2001). This requires a cut-off at the high- 
energy end of the distribution of the electrons, and adaptations have been made for such 
cases (Bhattacharya 2001; Panaitescu & Kumar 2001). The details of these adaptations, 
however, are still under debate, since the evolution of the high-energy cutoff is not well 
constrained. Since the lowest values for p we find are ~ 2. we do not take any adaptations 
for this effect into account. 

For several bursts the underlying parameter set has often been measured independently. 
Determinations for sets of GRBs have been made by e.g. Wijers & Galama (1999); Panaitescu 
& Kumar (2001); Yost et al. (2003); Gendre & Boer (2005); Granot et al. (2006). Such studies 
are generally limited to part of the parameter set, since they often use only one waveband, 
and therefore lack the possibility of finding the characteristic break frequencies in the broad- 
band spectrum. Unfortunately, there are few events where measurements in all wavebands 
are available, most notably GRB 970508 (e.g. Wijers & Galama 1999; Panaitescu & Kumar 
2002; Yost et al. 2003, Van der Horst & Wijers in preparation) and GRB 030329 (e.g. Berger 
et al. 2003b; Willingale et al. 2004; Frail et al. 2005). 

Here, we fit the broad-band spectral energy distributions (SEDs, from nIR through X- 
ray) of a subsample of the BeppoSAX sample of GRB afterglows. We constrain a subset of the 
blastwave parameters, namely the index of the power law energy distribution of electrons, p, 
the density profile parameter of the circumburst medium, k, and for some bursts the position 
of the cooling break, u c . Because of the paucity of radio data for most bursts in our sample, 
we have not included these, ensuring a more homogeneous approach. 

We make use of simultaneous fits in count space to obtain the most accurate possible 
measurements. In Paper I we provide details of the observations, data reduction and fitting 
method, summarized here in Sections 2 and 3 together with description of the model used. 
In Section 4 we present the results of our p- and A;- value and cooling break analysis, both 
for the sample and for individual sources. We compare these results to those of previous 
studies of this kind in Section 5, and perform statistical modeling of the p- value dataset. 
Here we discuss our findings in the context of the fireball model and long GRB progenitors. 
We conclude by summarizing our results in Section 6. 
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2. Observations and Spectral Fitting Method 

X-ray observations were made with the narrow field instruments on-board BeppoSAX 
(Paper I, Table 2), and here we have combined data from the MECS2 and MECS3 instru- 
ments (except in the case of GRB 970228, where we use the MECS3 instrument only). 

Optical and nIR photometric points were taken from the literature (Paper I, Table 3) 
and from our own observations of GRB 990510 (Curran et al. in preparation). Temporal 
decay slopes were again taken from the literature (see Paper I, Table 1): the temporal optical 
slopes from Zeh et al. (2006) and the X-ray temporal slopes from Gendre & Boer (2005); 
in ’t Zand et al. (1998); Nicastro et al. (1999). We have transformed the optical and nIR 
photometry to the log of the midpoint of the combined X-ray observations. We avoid using 
data taken before 0.1 days after the GRB, hence we assume no complex and flaring behavior 
occurs in the lightcurves and we restrict ourselves to the slow cooling regime. All data are 
transformed to count space for fitting purposes, in order that no model need be assumed 
a priori for the X-ray spectrum, and fitted within the ISIS spectral fitting package (Houck 
& Denicola 2000). Models consist of either a single or a broken power law, to allow for a 
possible cooling break in between the optical and X-rays. In the broken power law model 
we fix the difference in spectral slope to 0.5, as expected in the case of such a cooling break. 
Both Galactic and intrinsic absorption are components in the models, allowing for either 
Milky Way, Large Magellanic Cloud (LMC) or Small Magellanic Cloud (SMC) extinction 
laws for the GRB host galaxy, at either Solar (Z 0 ), LMC (Z= 1/3 Z 0 ) or SMC [Z— 1/8 
Z 0 ) metallicity. All details of the observations, reduction and analysis are given in Paper I, 
together with the results of the power law plus extinction fits. 


3. Theoretical Modelling 

We assume that the ambient medium density as a function of radius can be described as 
a power law with index i.e. n(r) oc r~ k . so that a homogeneous medium is given by k — 0 
and a stellar wind environment by k = 2 — the two most likely scenarios. For the purpose 
of looking at optical and X-ray data at ~ 1 day, we need to derive the time dependency of 
the peak flux F^ max and the cooling frequency v c as a function of k. : and the peak frequency 
v m (which has no dependence on k) , assuming the afterglows are in the slow cooling regime. 
After the jet break these parameters have no dependence on k ; in this region we know we 
are dealing with a jet geometry hence we label this case ‘Jet’ or j. See Table 1. 

If one assumes that the flux is a power law in frequency and time with /3 (or P) the 
spectral slope and a the temporal slope, using the conventions F u oc v ~H ~ a oc v ~ a , 
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with power law photon index T = (3 4-1, one can derive these slopes as a function of k and 
the power-law index p of the electron energy distribution. These values for a , (3 and T are 
given in Table 2 for two different situations: the observing frequency in between u m and u c , 
and the observing frequency above both frequencies. Also shown in this table are the closure 
relations between a and (3. 

One can invert all these relations to obtain p from a and /?, and even determine k from 
these observables: 

= 4(3(3 -2a) = 4[3(r — 1) — 2a] 

3/? — 2a — 1 3r - 2a - 4 1 J 

From Table 2 it is clear that when the observing frequency is higher than u m and u c , 
both a and (3 only depend on p, and do not depend on k. In the situation where the observing 
frequency is situated in between u m and u c , the spectral slope only depends on p, but the 
temporal slope depends on both p and k. So the structure of the ambient medium can only 
be determined in the latter case (v m < u < u c ), although having more accurate information 
on p from the situation with v m < v c < v is useful to get a better handle on k. 


4. Results 

The results of fits to the SEDs for all GRBs in the sample are given in Paper I in both 
Table 4 and Figure 2. For derivation of the blastwave parameter p we adopt the best-fitting 
models as listed in Paper I and in particular cases additional models were included. SMC- 
like absorption was the preferred extinction model in all cases except for GRB 000926 where 
LMC-like absorption is the best-fitting model. We calculate the values for p and k for two 
cases: the cooling frequency in between the optical and X-rays, and the cooling frequency 
above both. We have checked whether the cooling frequency could lie below the optical 
band using the relations of the fireball model, but this turns out not to be the case for these 
GRBs. The cooling frequency, v c , is obtained directly from the SED fits for 3 bursts: GRBs 
990123, 990510 and 010222, with u c = 1.31J5 x 10 17 , 4.3 ± 0.5 x 10 15 , 4.1^ 2 x 10 15 Hz 
at 1.245, 1.067 and 1.511 days since burst respectively. Applying the fireball model we find 
that another two sources, GRBs 980329 and 980703, may require a break within their SEDs, 
at ~ 2.6 x 10 17 and ux — 3 x 10 17 respectively. 

The resulting values for p can be found in Table 3 and for k in Table 4. All errors are 
given at the 90% confidence level for one interesting parameter, unless otherwise stated. a 1 
and a 2 refer to the pre- and post-break optical lightcurve slopes given in Paper I; we allowed 
for the possibility that these breaks are not jet breaks by considering that a 2 'is both pre- 
columns 6 and 7) and post- (column 8) jet break. This has also been done for the X-ray 
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Table 1: The temporal power law indices of the peak flux F Ujmax , the cooling frequency v c 
and the peak frequency u m as a function of k for pre- (columns 2-4) and post- jet-break 
(column 5). 
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Table 2: The temporal and spectral slopes of the flux, a and (3 (or T, where T = (3 + 1) 
respectively, the closure relations between a and (3 (or T), and p as a function of a, [3 and T. 
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temporal slopes. 


4.1. Results: individual sources 

4.1.1. GRB 970228 

We find p = 2.44^;^ and k = 1.731^695 w ^h u x < v c at 0.52 days. A cooling break 
between optical and X-ray bands is not required at the time of our SED (0.52 days) according 
to the F-test: the F-test probability, the probability that the result is obtained by chance, 
is 2.1 x 10~ 2 , which is quite high; so adding one free extra parameter is not a significant 
improvement. Using the best fitting model of a single power law plus SMC-like extinction, 
in the regime ux < v c , we find the data can be fit by both a homogeneous and a wind-like 
circumburst medium; the value of k is best constrained by the optical temporal slope. A 
break in the 77-band lightcurve has been reported at ^ 6 days (Galama et al. 1997), well 
after the time range covered by our SED. 


4.1.2. GRB 970508 

We find p = 2.56^oie an( ^ ^ = 0 . 49 +o; 67 i with < v c at 1.68 days. Using the single 
power law with LMC extinction (rather than SMC to obtain the more conservative errors 
on the spectral slope) and optical to X-ray offset free, in the regime v x < z/ c , we find the 
data are best fit by a homogeneous medium; the value k is best constrained by the optical 
temporal slope. Broadband modeling by Van der Horst & Wijers (in preparation) with k as 
a free parameter gives very tight constraints on k: a value of 0.02 ± 0.03 is derived, i.e. a 
homogeneous medium. 

Previous works put the cooling break at optical frequencies, v c — 1.6 x 10 14 Hz, at 12.1 
days since burst, between the B and V bands (e.g. Wijers & Galama 1999). We, however, 
find that the cooling break is likely to lie above the X-rays at 1.68 days. We note in this 
context the uncertain extrapolation of the optical data used in the SED, owing to an irregular 
shaped lightcurve at early times, which we have attempted to account for in allowing the 
optical to X-ray offset to go free. This is a particularly difficult case given that the X-ray 
data cover the time period immediately following an optical flare when the optical lightcurve 
appears to have flattened before breaking to a typical and well defined power law. 
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4-1.3. GRB 971214 

We find p — 2.20 ± 0.06 and k = 2.17lo!35 ? with v x < v c at 1.36 days. Using the 
best-fitting model of a power law plus SMC extinction, in the regime vx < v c , the data are 
best fit by a wind medium: from the optical temporal slope we find k — 2.1 7^35, whilst 
from the X-ray temporal slope k = 2.33tolt- A broken power law provides no acceptable 
fits, i.e. the F-test probability is high (4.1 x 10~ 2 ). A spectral break is claimed for this burst 
in the IR (~1 micron) at 0.58 days (v c ~ 3 x 10 14 Hz, Ramaprakash et al. 1998). This is 
not the cooling break, but the peak of the SED moving to lower frequencies, so there is no 
conflict with our results for v c . 

4 .I. 4 . GRB 980329 

We fin d p = 2.50^02 an d & = — 4.891 25 4 q 8 , with u c < vx at 1.15 days. The spectral 
fit obtained with a single power law plus SMC extinction is inconsistent with the optical 
temporal slopes. A spectral break in the power law is not favoured by the F-test (probabil- 
ity of 7.2 x 10 -2 ), but this spectral break model provides agreement between the spectral 
parameters and the optical and X-ray temporal slopes. In this regime, v c < ux, k cannot be 
well constrained though the centroid of the fit to the optical data is that of a homogeneous 
medium. We note that when omitting the I band point from the SED. which may be over- 
estimated (see Yost et al. 2002, and Paper I), our results do not change. In the absence of 
a spectroscopic redshift determination for GRB 980329 we adopt the photometric redshift of 
2 = 3.6 (Jaunsen et al. 2003), hence all results should be treated with caution. 


4.1.5. GRB 980519 

We find p = 2.96^;Qg and k — 0.23U,; 22 , with vx < v c at 0.93 days. Using the power 
law plus SMC extinction model for the SED, in the regime isx < v c , the optical data are 
best fit by a homogeneous medium, and the X-ray temporal slopes can be fitted by both a 
homogeneous and a wind medium; k is therefore best constrained by the optical temporal 
slope. In contrast, Chevalier & Li (1999) found that the radio emission of the afterglow of 
GRB 980519 measured between 7.2 hours and 63 days since the burst is consistent with an 
external wind instead of a homogeneous medium. Frail et al. (2000) note, however, that the 
interstellar scintillation present in the radio data does not allow to draw firm conclusions on 
this. 

The optical temporal break at 0.48 ± 0.03 days (Zeh et al. 2006) cannot be explained 
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by passage through the optical bands of z/ c , since the derived value for p from the temporal 
slope is too high in that case (p = 3.69 ± 0.06) compared to the p-value from the joint 
spectral slope. It also cannot be explained by a jet break, since p is too low in that case 
(p — 2.27 ± 0.05). It appears that the fireball model is a good explanation for the first 
temporal slope and the spectrum used here, but the post-break optical slope has either been 
incorrectly measured or we do not yet have the correct model for this afterglow. 

We caution that all the results are based on a redshift estimate (see Paper I) equal to 
the mean of the sample spectroscopic redshifts of z — 1.58, but an independent estimate 
of the redshift using the X-ray lightcurve (B. Gendre, 2006) also provides a value of the 
same order (z = 1.4 ± 0.2). This afterglow showed a very steep temporal decay compared 
with other GRBs (Halpern et al. 1999). This is somewhat reminiscent of the very early-time 
decays of many Swift bursts which occur at <500 s after trigger (Nousek et al. 2006) and 
likely have a significant prompt emission component. 


4.1.6. GRB 980703 

We find p — 2.74^48 and & = with u c < z/* at 1.33 days. The spectrum is 

best fit with a single power law plus SMC extinction. This is, however, inconsistent with the 
X-ray temporal slopes. The broken power law model, though not favoured according to the 
F-test (probability of 7 x 10~ 3 ), is necessary to obtain agreement with the X-ray temporal 
slopes. We use the broken power law plus SMC extinction model, in the regime v c < vx- 
The nature of the optical temporal break at 1.35 ± 0.94 days cannot be determined because 
of large uncertainties in the optical temporal slopes, which are also the reason why k cannot 
be constrained. 

Two publications have postulated a position for the cooling break in past studies. 
Vreeswijk et al. (1999) propose v Q < vx < v c. at 1-2 days after the burst, Bloom et al. 
(1998) propose v Q < u c < vx at 5.3 days, and our SED study at 1.33 days, when compared 
with the optical temporal slope, is inconclusive since both vx < v c and vx > v c can be 
accommodated. It may be that the cooling break has moved into our observed bands during 
accumulation of the X-ray spectrum (possibly indicated by the inability of the fireball model 
to fit the data when a single power law is assumed for the spectrum). If we require consis- 
tency with these previous results, the cooling break must be moving to lower energies and 
lies approximately at X-ray frequencies in our data. This would mean that the circumburst 
medium is homogeneous, since the u c is expected to move as £~ 1</2 in this case, while z/ c will 
increase in time as t 1 ^ 2 in the wind case. 
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The host galaxy of GRB 980703 appears to have a high and possibly variable optical 
extinction along the line of sight to the GRB (see Paper I). The different (and formally 
inconsistent) values of Ay may be due to different methods for measuring the extinction, 
probing of different regions of the host galaxy, or may indicate that the environment in which 
the burst occurred is changing with time. We have used the Vreeswijk et al. (1999) optical 
data and scaled it from 1.2 days to 1.33 days after trigger. Combining the optical and X-ray 
data when fitting provides us with a different estimate for the extinction than was obtained 
by Vreeswijk et al. (1999) for the optical alone. Any change in optical extinction will have 
an effect upon the measured spectral slope and hence the derived value of p. 


4.1.7. GRB 990123 

We find p = 1 -99^07 an d k = 2.00loJi! w *th at 1.25 days. Using the best 

fitting model of a broken power law plus SMC extinction, in the regime < v X -, the optical 
data are best fit by a wind medium. The optical temporal break at 2.06 ± 0.83 days is 
marginally consistent with a jet-break: the p- value derived from the post-break temporal 
slope is consistent with the one derived from the pre-break optical temporal slope, but 
inconsistent with the spectral slope. The uncertainties in the X-ray temporal slope are too 
large to determine the phase of blastwave evolution (i.e. before or after jet-break). 

For GRB 990123 the temporal slope difference between optical and X-ray of 0.25 also 
agrees with the spectral analysis where a broken power law model is the best fit, indicating 
a cooling break between the optical and X-ray bands at 1.25 days post-burst. The value we 
derive for p is consistent with that derived from the X-ray spectrum alone of p — 2.0 ±0.1 
(Stratta et al. 2004) and lower than a previous estimate via broad-band modelling of p — 
2.28 ± 0.05 (Panaitescu & Kumar 2002) (we note that the latter uncertainty is 1 a and not 
the 90% error used in the rest of our paper). 


4.1.8. GRB 990510 

We find p = 2.06^02 an d & = O.SO^go, with u c < at 1.07 days. The best fitting 
model to the SED is clearly a broken power law with negligible extinction, in the regime 
u c < v x . There is considerable improvement in the y 2 when allowing for this break in the 
power law, noted by previous authors, which we find is located at 0.016 - 0.020 keV at ~1.07 
days since burst (consistent with the value given by Pian et al. 2001). The slope change is as 
expected for a cooling break in the slow cooling regime when leaving both power law slopes 
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free. 

The optical data are best fit by the values for p and k mentioned above. In this case, 
however, the X-ray temporal slope is not consistent with the spectral slope nor the optical 
temporal slope at the 90% (~ 1.6or) level, but does agree within 3 a. The optical temporal 
break at 1.31 ± 0.07 days is marginally consistent with a jet-break. The derived value of p 
is consistent with the value derived from the BeppoSAX X-ray spectrum alone of p ~ 2.1, 
by Kuulkers et al. (2000), and rules out the value of p = 2.6 ± 0.2 also derived from the 
BeppoSAX X-ray spectrum by Stratta et al. (2004); that same X-ray data is used here, but 
is combined with nIR and optical data to obtain our limits on p. 


4.1.9. GRB 000926 

We find p = 2.54^o;Qg and k = 2.16±oio! with zqc < z'c at 2.23 days. Using the power 
law plus LMC extinction model, i.e. in the regime r'x < z'c, the optical data are best fit by 
a wind-like medium. Large uncertainties in the X-ray temporal slope prevent determination 
of the circumburst medium structure and blastwave evolution phase; the optical temporal 
break at 2.10 ± 0.15 days is consistent with a jet-break. 


4.1.10. GRB 010222 

We do not find a consistent solution for this afterglow taking 90% uncertainties, but we 
do find one taking 3 a uncertainties: p = 2.041q^o and k = 2.28^o!293 with u c < vx at 1.51 
days. Adopting the single power law model with LMC extinction, in the regime vx < v c , 
the optical slopes are not consistent with the spectral slope nor the X-ray temporal slopes, 
both at 90% and 3(7 level. Using a broken power law plus SMC extinction, in the regime 
< vx, for which the F-test indicates a marginal improvement (probability of 1 x 10~ 4 ), 
the X-ray temporal slope is not consistent with the spectral slope nor the optical temporal 
slope at 90% level, but they are consistent at 3o\ In the first case we obtain p = 2. 721^ 
and k — O.SO^oJ®; in the latter case p = 2.04loao and k = — 2.25I43 04 , which is derived 
from the pre-break optical temporal slope, or k = 2. 281^29, derived from the post-break 
optical temporal slope. Since the temporal break happens quite early, 0.64 ± 0.09 days after 
the burst (Zeh et al. 2006), and the post-break optical slopes are inconsistent with a jet- 
interpretation, the early temporal slope is probably influenced by late-time energy injection 
and a medium with k = 2.281 q 29 th e correct interpretation. 
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Panaitescu & Kumar (2002) find the cooling break to lie at optical wavelengths or 
longer, in agreement with our results. However, they derive a low value for p of 1.35, 
and find significant reddening of the optical spectrum which they say explains the second 
steepening observed in the optical after 6 days by Fruchter et al. (2001). Bjornsson et al. 
(2002) argue that the unusually slow decay of this afterglow and positive detection of linear 
polarization can be explained by a jet model with continuous energy injection. Such slow 
decays are seen in the ‘canonical’ Swift X-ray lightcurves (e.g. Nousek et al. 2006), termed 
the plateau phase, and typically begin a few hundred seconds after the GRB trigger. The 
electron energy distribution then has p — 2.49 ± 0.05, which is inconsistent with our result 
for p from the spectral fits. 


5. Discussion 

The parameters that can be derived from broadband modeling of afterglow lightcurves 
describe the micro- and macrophysics of the relativistic jet and its surrounding medium. To 
obtain the full set of parameters the spectral energy distribution has to be covered from X- 
ray to radio wavelengths. Two of these parameters can be deduced from just the spectral and 
temporal slopes in the optical and X-rays. i.e. the electron energy distribution index p and 
the circumburst medium profile parameter k. These two parameters have been determined 
in this paper for a selection of 10 GRBs, for which the final results are shown in Table 5 and 
Figure 1. For completeness the values for v c are also listed in Table 5; half of the GRBs in 
this sample have a value v c < Px- 


5.1. The distribution of p 

Some theoretical studies of particle acceleration by ultra-relativistic shocks indicate that 
there is a nearly universal value of p of 2.2 — 2.3 (e.g. Achterberg et al. 2001), while other 
studies indicate that there is a large range of possible values for p of 1.8 — 4 (e.g. Baring 
2004). From the results presented in this paper and from broadband modeling of individual 
bursts, quite a large range of values for p have been found, which could indicate that there is 
a large intrinsic scatter in the value of p. Here we test the null-hypothesis that the observed 
distribution of p can be obtained from a parent distribution with a single central value of p, 
by performing a statistical log-likelihood analysis on the obtained values of p. 

We generate synthetic datasets for p from the 10 bursts in our sample by taking random 
numbers from the normal distributions that are described by the most likely value of p and 



Table 3. Values for p. We calculate the results for the cases k — 0 and k = 2. Bold type highlights consistent results 
(at 90% level). In cases where the best-fitting spectral model to the SED (from Paper I) is inconsistent with the 
model fits, we show the results for this original best- fitting model in italics. 
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Table 4: Values for k. Bold type highlights consistent results (at 90% level). In cases where 
the best-fitting spectral model to the SED (from Paper I) is inconsistent with the model fits, 
we show the results for this original best-fitting model in italics. 
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Fig. 1. — Derived values of p (top panel) and k (lower panel) for each individual afterglow 
(see Table 5): the horizontal axes represent the GRBs in date order left to right and errors 
are 90% confidence. In the top panel the dotted line indicates the most likely value of 
p = 2.035; in the lower panel the dotted lines indicate k = 0 (homogeneous medium) and 
k — 2 (stellar wind). 
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Table 5: Final results for p, k and v c for all ten bursts in the sample. 


“Optical data extrapolation is uncertain 

^Broken power law better than single power law in SED. 

c This solution is consistent with all measurements except the post-break optical temporal slope. 

^ Broken power law better than single power law in SED. Large uncertainties in the optical temporal slopes. 
e X-ray temporal slope only consistent at 3 a level. 

^ X-ray temporal slopes have large uncertainties. 

9 Break in the optical lightcurve at 0.6 days is not a jet-break. X-ray temporal slope only consistent at 30 - 
level. 
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the uncertainties in the 10 values of p. We take the asymmetry of the errors into account by 
demanding that the two halves of the distribution are connected at the p value and that the 
total probability is equal to 1. This results in normal distribution functions which are given 

by 


f(x-,p,a 1 ,a 2 )dx = 


v/2 


v^Oi + o- 2 ) 


dx 


e -(x-p) 2 /2al 

-{x-p) 2 /2crl 


(x < p ) 
{x > p ) 


( 2 ) 


with <ti and a 2 the lower and upper 1 a uncertainties in p respectively. To convert the 90 % 
confidence limits in table 5 to 1 a uncertainties, we divided those by a factor of 1.6. 


The log-likelihood for asymmetric normal distributions is given by 


-2 log A = -2 log (jj/i ) = ^ (-2 log £) 
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log 


<?l,i + a 2,i 
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in which in the last term a lti is used for x < Pi and <j 2 ,i for x > pi- We determine the 
minimum of the log-likelihood for our 10 bursts and find the most likely value of p — 2.035 
with a log- likelihood of 613.6. We generate 10 4 synthetic datasets from the most likely p- 
value and the 1 a uncertainties in the 10 values of p. We then take 10 values of p coming 
from these synthetic datasets (one random number from each of the 10 datasets), calculate 
the most likely value of p and the accompanying log-likelihood. The resulting distribution of 
log-likelihood values is plotted in figure 2, together with the minimal log-likelihood coming 
from the real data. We find that in 99.92 % of the cases the log-likelihood of the synthetic 
data is smaller than the one coming from the real data. This leads to the conclusion that 
the hypothesis, that the distribution of p from our sample can be obtained from a parent 
distribution with a single central value of p, is rejected at the 3 cr-level. 


From Table 5 and Figure 1 it can be seen that the log-likelihood is dominated by the 
value for p of GRB 990123, and after that by GRB 980519. We have performed the same 
analysis after leaving out one of the values for p and find that in 8 of the 10 cases the 
hypothesis is rejected at a level of > 3<r, including the omission of GRB 010222 for which 
the closure relations were not all well satisfied. Only in the cases of leaving out GRB 990123 
or GRB 980519 we find that the hypothesis is rejected at a ~ 2 a- level. We note that taking 
only the p - values of GRB 990123 and GRB 980519, we find that the hypothesis is completely 
rejected. 

If the microphysics of each burst is essentially the same, one would expect p to be 
constant. The value p = 2.2 has been widely quoted as a typical number and the deviations 
from this interpreted as due to the external environment or further energy injections from 
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log-likelihood 

Fig. 2. — The log-likelihood distribution coming from the 10 4 synthetic datasets generated 
from one single value of p (solid line); the dashed vertical line indicates the log-likelihood for 
the real data. In 99.92 % of the cases the log-likelihood of the synthetic data is smaller than 
the one coming from our measured sample. 



-19- 


the source (e.g. Berger et al. 2003a). Meszaros et al. (1998) showed that for a single value 
of p , variations in the jet energy per solid angle, i.e. structured jets, could lead to a range 
of lightcurve decays; which can also be produced if a set of intrinsically similar, structured 
jets are viewed at various off-axis angles (Rossi et al. 2002). The study of such structure in 
the jet and viewing angle dependence is, however, beyond the scope of this work. 

Johannesson et al. (2006) claim that from broken power law fits the value of p is overes- 
timated from the pre-break lightcurve slope while being underestimated from the post-break 
slope. The results presented in that paper were however only for a homogeneous ambient 
medium; they claim that in a wind medium the p value from the pre-break lightcurve slope 
does not show a systematic deviation. 

More recently, breaks in lightcurves, both optical and X-ray, which would generally have 
been attributed to jet-breaks, have been found to be achromatic, in sharp contrast to the 
picture of a jet-break: the sideways expansion is not expected to be frequency dependent. 
As such, several previously claimed jet-breaks may have to be revised; unfortunately, there 
is no relevant data to confirm or reject these claimed jet-breaks, since these have been found 
only in one waveband (mostly the optical, owing to the lack of dense monitoring in X-rays 
in the pr e- Swift era). As such, the mention of jet-breaks in this paper has been taken at face 
value, but with the caveat mentioned here. 

If we compare the values for p calculated here with those calculated from the X-ray 
spectra alone (Stratta et al. 2004), we find they agree within the 90% errors except for the 
bursts GRB 970508 and GRB 990510. With this method we reduce the average 90% error 
on p from ±0.58 (Stratta et al. 2004) to ^o' 20 ) an d the values obtained are more accurate 
given that consistency between the nIR, optical, UV and X-ray is required. 

Panaitescu & Kumar (2002) calculated jet parameters for a sample of 10 GRBs including 
several studied here, using broadband observations including radio data when possible. In 
their model, based also on the fireball model, they assume uniform jets (no structure across 
the jet) and the the energy parameters e e and e# are constant, and finally they assume the 
observer is located on the jet symmetry axis. Our spectrally-derived values are consistent 
with theirs at the 3cr level or better for GRBs 970508, 980519, 990510 and 000926. There is 
no agreement for GRB 990123. 

They found a spread in p values as do we, but with half the values lying below 2, and 
a mean value of p = 1-92^26 (2cr). 

Chevalier & Li (2000) carried out a similar study of broadband afterglow data, and their 
estimates for p are in agreement with ours for the sources 970228, 970508, 980519 and 990510, 
and disagree for 990123 (they do not quote errors per source but estimate errors to be ~0.1). 
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They conclude that their sample shows a range in the values of p which is not consistent 
with a single value. 


5.2. The circumburst medium 

The profile of the circumburst medium is a particularly important parameter in studying 
the progenitors of gamma-ray bursts. In the case of long-soft bursts the progenitor is a 
massive star that is expected to have had a massive stellar wind in earlier phases of its 
evolution, and one would expect to see a signature of that in the afterglow lightcurves. 
Evidence for a stellar wind in the form of fast outflowing absorption lines within restframe 
UV spectra has been seen in a small number of cases, the best example being GRB 021004 
(e.g. Schaefer et al. 2003; Starling et al. 2005). This does not mean, however, that a density 
profile with n oc r~ 2 is always expected, since this assumes a constant mass-loss rate and a 
constant wind velocity. Changes in mass-loss rate and also interactions of the wind with the 
interstellar medium can alter the ambient medium profile. 

In the first broadband modeling attempts the ambient medium was assumed to be the 
homogeneous interstellar medium, which was consistent with the derived particle densities. 
However, since the progenitors of at least a fraction of all GRBs are now known to be massive 
stars and the blastwave is situated at ~ 10 17 cm at ~ 1 day after the burst, a massive stellar 
wind profile is expected. Nonetheless, the medium that the blastwave is probing could still 
be homogeneous, because of the emergence of a reverse shock in the wind when the wind 
meets the interstellar medium (see e.g. Wijers 2001; Ramirez- Ruiz et al. 2005). This shocked 
wind turns out to become homogeneous and, for typical physical parameters derived from 
afterglow modeling, the blastwave encounters the transition from the stellar wind to this 
homogeneous shocked medium at ~ 1 day (see e.g. Pe’er & Wijers 2006). The actual time of 
the transition, which would be detectable in the afterglow lightcurves, depends for instance 
on the mass-loss rate and the density of the interstellar medium, which are both not really 
well constrained for most GRBs. In our sample we do not see this kind of transition in 
the optical lightcurves in which there is a break. It has been claimed for some bursts, for 
instance GRB 030329 (Pe’er & Wijers 2006), that this transition is observed, but the bumpy 
lightcurves cause confusion. 

Another way to obtain a constant density from a massive stellar wind is in the region 
after the wind termination shock. The distance to the termination shock can be very large, 
but recent observations of two Wolf-Rayet binaries has suggested that this distance could be 
several times smaller if the wind is asymmetric. Eldridge (2007) shows that wind asymmetry 
probably exists in two systems, that can be caused for example by rotation which is expected 
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for GRB progenitors in the framework of the collapsar model in order to retain enough 
angular momentum. If the asymmetry exists for the entire stellar lifetime, then a closer 
termination shock and asymmetric supernova may be expected, increasing the chances of 
observing an afterglow traversing a constant density medium. 

In our sample of 10 GRBs there are four sources that are consistent with an r~ 2 wind 
medium, with relatively small uncertainties, namely GRB 971214, GRB 990123, GRB 000926 
and GRB 010222. There are three GRBs which are not consistent at a 90 % confidence level 
with a wind medium, GRB 970508, GRB 980519, and GRB 990510, although for the first, 
one caution is warranted with the interpretation of the lightcurves; and for the other four 
bursts the uncertainties are too large to constrain the ambient medium profile. We contrast 
our findings with Panaitescu & Kumar (2002), who, in broadband fits to the data of 10 
bursts, found that a wind-like medium was preferred over a homogeneous medium in only 
one case: GRB 970508. For this particular burst our analysis provides a value of p which is 
an equally good fit to wind or uniform medium predictions for u < u c from the spectra and 
lightcurves, but the closure relations are obeyed ( at the 2<r level for both cases) if v > v c . 
They find circumburst densities of order 0.1-100 cm -3 for most sources, which they say 
demonstrates the surrounding medium does not have, in general, the r~ 2 profile expected 
for the unperturbed wind of a massive GRB progenitor. 

The association of long-soft gamma-ray bursts with Ib/c supernovae was first suggested 
for GRB 980425 / SN 1998bw by Galama et al. (1998), and confirmed by the discovery of 
GRB 030329 / SN 2003dh (e.g. Hjorth et al. 2003). Therefore, it is useful to compare the 
circumburst medium characteristics derived from gamma-ray burst afterglows and from radio 
observations of supernovae, which also trace the density profile of the surroundings of these 
massive stars. From radio supernovae r~ 2 density profiles have been found, but also in some 
cases significantly flatter behaviour of ~ r -1 - 5 [ n SN 1993J and SN 1979C (for a review on 
radio supernovae see Weiler et al. 2004, and references therein for individual supernovae). 
In the latter case a transition of r ~ 2 to r~ 1A was even observed. This flatter density profile 
can be attributed to changes in the mass-loss rate of the massive star in some phases of 
its evolution. The three sources in our sample of GRBs with a density profile flatter than 
r~ 2 are possible examples of the relativistic blastwave ploughing its way through a region 
of the circumburst medium which is affected by changing mass-loss rates. So although in 
Table 5 we have described them as GRBs with a homogeneous ambient medium, this is not 
necessarily the case. Especially for GRB 970508 and GRB 980519 this is a possibility, but 
the uncertainties on k are too large to claim whether it is a homogeneous or flattened wind 
medium. GRB 990510 has smaller uncertainties and seems to be more consistent with a 
homogeneous medium, especially since the upper limit on k is ~ 1.0, which is much flatter 
than what is observed in radio supernovae. 
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The Swift satellite now provides us with substantially greater coverage of a large number 
of X-ray afterglows (~100 per year) and often with high quality data from which to measure 
the spectral and temporal slopes. However, few of these also have substantial optical follow- 
up. The combination of X-ray and optical data helps determine the position of the cooling 
break and obtain accurate spectral slopes which provide the value of p. For the derivation of 
k in this study, we have found the optical temporal data most constraining. For this reason, 
and for the confirmation of achromatic jet breaks it is essential that such optical data be 
obtained for as many well sampled Swift X-ray afterglows as possible. 


6. Conclusions 

We have measured the distribution of the the injected electron energy distribution index, 
p, and the density profile of the circumburst medium, n(r) oc r -fc , from simultaneous spectral 
fits to their X-ray, optical and nIR afterglow data of 10 BeppoSAX GRB afterglows. 

We find a most likely value of p for the sample as a whole of 2.035. A statistical analysis 
demonstrates that in fact the distribution of p values in this sample is inconsistent with a 
single value for p at the 3<r certainty, which is at variance with many theoretical studies. 

We measure the distribution of the local density parameter k, generally only assumed to 
be 0 or 2, and we find that the majority of GRBs for which we can constrain k are consistent 
with a wind-like circumburst medium. One source is clearly, i.e. > 3cr certainty, inconsistent 
with this picture and fits instead a homogeneous medium. These results are consistent with 
the expectations of at least a subset of GRBs originating from massive stars, which have 
a lot of mass-loss in the form of a surrounding stellar wind. We discuss the possibility of 
values of k < 2 within the stellar wind framework. 

The method presented here provides a way to study the distribution of blastwave pa- 
rameters for a sample of GRBs, and allows estimates to be derived when insufficient data are 
available for a full time-dependent solution. In the current Swift-era the method is equally 
well applicable, although one would have to ensure that the data are in the afterglow domain, 
i.e. not contaminated by prompt emission or late-time energy-injection. A decent sampling 
of the optical SED and lightcurve, more difficult with the average fainter Swift afterglow 
sample, is crucial to constrain the temporal slopes and cooling break frequency. 
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